Non-substitutional single-atom defects in the Gei_^Sn^ alloy 
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Gei-i^Siii, alloys have proved difficult to form at large x, contrary to what happens with other 
group IV semiconductor combinations. However, at low x they are typical examples of well-behaved 
substitutional compounds, which is desirable for harnessing the electronic properties of narrow 
band semiconductors. In this paper, we propose the appearance of another kind of single-site defect 
(/9 — Sn), consisting of a single Sn atom in the center of a Ge divacancy, that may account for these 
facts. Accordingly, we examine the electronic and structural properties of these alloys by performing 
extensive numerical ab-initio calculations around local defects. The results show that the environ- 
ment of the (3 defect relaxes towards a cubic octahedral configuration, facilitating the nucleation 
of metallic white tin and its segregation, as found in amorphous samples. Using the information 
stemming from these local defect calculations, we built a simple statistical model to investigate 
at which concentration these (3 defects can be formed in thermal equilibrium. These results agree 
remarkably well with experimental findings, concerning the critical concentration above which the 
homogeneous alloys cannot be formed at room temperature. Our model also predicts the observed 
fact that at lower temperature the critical concentration increases. We also performed single site 
effective-field calculations of the electronic structure, which further support our hypothesis. 

PACS numbers: 71.20.Nr,71.55.Ak,71.15.-m 
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I. INTRODUCTION 



The search for direct energy-gap materials based on 
group IV semiconductors is not newii^i^i^, but has been 
hindered for a long time by sample preparation problems. 
The importance of such compounds due to their poten- 
tial for technological applications was recognized early. 
A low-energy direct gap would enable the use for in- 
frared applications, long-wavelength photodetector sand 
emiters, LED and infrared laser applications, i.e. in gen- 
eral allow the use for optoelectronic devices. Combined 
with the use of group IV elements, this would be ideal 
for compatibility and integration with the Si-based tech- 
nology. 

While Ge is an indirect-gap semiconductor (0.41 eV 
gap, fundamental energy band-gap is 0.76 eV at the di- 
agonal L symmetry point of the fee diamond Brillouin 
zone)^, gray tin (a — Sn) is a semimetal with the same 
crystal structure (no band gap, at the zone center its 
valence and conduction bands overlap about 0.42 eV)^. 
It was recognized that Gei-^^Snj^ alloys should provide 
tunable gap materials, controlled by concentrationj^i^i^ 
and a route for obtaining direct-gap group IV systems 
compatible with Si integrated circuits. The first exper- 
iments providing evidence for the indirect-to-direct gap 
transition date from 1997^ while only recently crystalline 
samples of the quality and stability required by device 
applications were preparedJ^iiiii^ In the last two years 
further important uses for these materials were proposed: 
they were studiedi^ in connection with their potential 
use for nanostrucured thermoelectric cooling devices, and 
the higher electron and hole mobilities which might be 
reached by alloying and strainingf^i Thus, they might al- 
low to overcome current integrated circuit limitations to 



develop higher speed Si microelectronics. 

The electronic structures of most semiconductor alloys 
are smooth functions of their composition, thus providing 
a versatile tool for device engineering. Alloys of elemen- 
tal semiconductors such as group IV elements Si and Ge, 
and alloys of III-V compounds as GaAs, AlAs,InAs and 
InP play key roles in microelectronics and optoelectron- 
ics. Sii-xGca; is nearly ideal as it has lattice constant 
and interband optical transition energies which are es- 
sentially linear in x. However, S'z-based photo-detectors 
e.g. cannot cover the optical communication wavelength 
windows (1310-1550 nm). Pure Ge barely reaches the 
1550nm range, but it doesn't grow well on Si substrates: 
this was solved by alloying Ge with Sn. Even 2% Sn 
achieves an order of magnitude increase of absorption at 
1550 nm, and the alloy grows with device quality on Si 
substrates J^iiiii^ii^ 

Binary alloys of group IV semiconductors are usually 
easy to prepare at any concentration, but this is not the 
case for the Gei^ajSux alloy. Homogeneous alloys have 
proved difficult to form at large x. This has considerably 
reduced its applicability in the fabrication of electronic 
devices. Other group IV alloys are formed by simple 
substitution of one component by the other, maintain- 
ing the overall tetrahedral environment of the diamond 
lattice. This is true even in amorphous alloys. 

At the low concentrations at which it has been possible 
to prepare homogeneous Gei-^Sua; alloySfiS it has been 
found that 5*71 occupies substitutional sites. This results 
in a smooth dependence of its electronic properties as a 
function of concentration, which is desirable to tailor key 
quantities as the band gap and the density of states at 
the Fermi level. However, at Sn concentrations higher 
than 20% this picture breaks down because Sn exhibits 
the tendency to segregate in the cubic f3 phase, which is 



2 



metallic. This spoils the semiconducting properties. 

The main point of interest of our present work is to try 
to understand theoretically the structural and dynam- 
ical mechanisms that trigger this peculiar segregation. 
Some authors^^ attribute this to the large difference in 
size between Ge and Sn atoms and to the large lattice 
mismatch. This cannot be the whole picture, because 
segregation has also been observed to begin at the same 
concentrations in amorphous alloysii. 

Our hypothesis is that the difficulty in the formation of 
homogeneous alloys can be explained by the nucleation 
of the metastable /3 — Sn phase ( "white tin" , which is a 
metal and cubic) in the perfect tetrahedral lattice of the 
alloy. The tetrahedral a phase, known as "gray tin" , is 
a semiconductor with a small direct gap. 

We shall investigate the mechanisms by which clus- 
ters of few Sn atoms can acquire a cubic symmetry in 
the tetrahedral Ge matrix. Accordingly, we assume that 
there exists a defect with a single Sn atom bonded to six 
neighbouring Ge atoms. This is possible due to the im- 
portant role of the d— orbitals in the electronic structure 
of both atoms. The existence of such a defect in amor- 
phous alloys has been confirmed by detailed Mossbauer 
experimentsii which show a signal corresponding to a 
Sn atom in an octahedral environment, besides the ex- 
pected signal of the substitutional tetrahedral a — Sn 
defect. Signatures of the presence of two sites for the 
Sn atoms were also found with ^^^Sn Mossbauer spec- 
troscopy in nanoscale mixtures of Ge/Sn}^ This point 
defect, that we shall call (3 ~ Sn defect from now on, can 
be envisaged as a single Sn atom at the center of a Ge 
divacancy. This defect (although rare) should produce a 
large negative (shrinking) elastic field around it, oppo- 
site to the usual positive (expansive) elastic field of the 
a — Sn defect, due to the larger size of Sn as compared 
to Ge. 

In the following section we will present our hypothesis 
for the incorporation of Sn in the Ge lattic in more detail, 
in Section III we will describe the different approaches we 
have used to tackle the problem sistematically as well of 
the main results obtained with each, providing a quite 
complete scenario and consistent explanation for the ex- 
perimental findings in the framework of the hypothesis 
we investigated. Finally we will present our conclusions, 
discussing how our results might be useful to solve prob- 
lems for sample preparation with device quality, as re- 
quired for the important technological applications en- 
visaged for these alloys, and give an outlook to further 
interesting research on these materials. 



II. HYPOTHESIS: MECHANISMS FOR Sn 
INCORPORATION IN THE Ge LATTICE 

The proposed mechanism goes as follows: small 
amounts of Sn incorporate very easily in isolated places 
in the Ge matrix as a defects. The strain caused by the 
size mismatch increases with the concentration of these 



defects (xa), causing a propitious environment for the 
formation of the (3 defect, since that would release the 
strain in the lattice. As the a defects get closer to the 
existing /3 defects they attract each other through their 
elastic fields, and if they merge, further elastic energy 
would be released and the small clusters of Sn atoms 
would relax to a cubic symmetry. Therefore, if difusion 
allows the migration of defects during the formation, the 
material would present a natural tendency to segregate 
the Sn clusters to its surface, in order to attain equilib- 
rium. 

All these ideas have to be tested either by numerical 
calculations or statistical models. First of all, one needs 
accurate ab-initio electronic calculations to estimate the 
scale of energies involved in the immediate environment 
around each defect. These calculations also provide the 
relaxed configurations of the defects, allowing one to test 
the cubic symmetry locally. Furthermore, by looking at 
small displacements within the clusters one can estimate 
the dependence of the electronic energy with the volume, 
which is a pressure directly related to the elastic field 
caused by the defects. All these quantities will allow to 
construct a model for the behavior of the macroscopic 
system, which will be needed to test if one can form /3 
defects in the alloy, and at which concentrations they 
would appear. 

The electronic properties of the alloy can be tested 
by performing effective- field calculations, like the virtual 
crystal approximation (VCA)^^ or the coherent potential 
approximation (CPA).^° It is well-known that these ap- 
proximations give excellent predictions in substitutional 
alloys, thus a failure to describe the experimental find- 
ings suggests a further indication of the presence of /3 
defects. In fact, we noticed that if one takes into ac- 
count the structural changes undergone as a function of 
Sn concentration, one can extend the range of validity of 
these approximations. 

In the following section we will present in detail these 
three approaches to the problem, and discuss the results 
in terms of the main hypothesis presented here. 



HI. NUMERICAL AND MEAN FIELD STUDIES 
A. Local defect electronic calculations 

A good starting point to examine the properties of 
local defects in an alloy is to perform ab initio elec- 
tronic calculations using the common techniques of den- 
sity functional theory. 

We have basically two different point defects to be 
studied: 1) a Sn atom sitting substitutionally in a per- 
fect Ge diamond lattice, and 2) a single Sn atom in a 
cage of six Ge atoms that surround a double vacancy in 
the perfect diamond lattice. This cage has a long axis 
along the (111) direction, but it is not clear if the pres- 
ence of the impurity inside the cage should modify this 
symmetry. 
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The electronic energies depend on the local environ- 
ment around the defect and the respective relaxation of 
the lattice, which involves changes in lattice parameter 
(or cell volume) and in the individual cell- atom positions. 
We have considered unit cells of different sizes around the 
point defect, in order to address these issues. 

On one hand, we analized to what degree the ener- 
gies and the lattice relaxation around the defect depend 
on the size of the cell, and estimated the elastic field 
caused by the defect. On the other hand, we studied 
more complex defect environments involving several Sn 
impurities. We have considered unit cells of 8-atoms, as 
the ones studied in previous works^i^, as well as some 
larger cells, with 16 and 64 atoms, designed to better 
assess the preferred configurations when multiple defects 
are present. 

We have determined the defect energies per atom, rel- 
ative to the pure Ge and Sn cohesive energy reference 
values, allowing full relaxation of the different cells stud- 
ied with respect to the Ge lattice (i.e. allowing changes 
of lattice parameter and of individual positions of the 
atoms in the cell) . Specifically, for a unit cell consisting 
of n Ge atoms and m Sn atoms, we define the energy of 
the defect (Ed) as the cohesive energy per atom: 

Ed — ^^Ccll(Gc„Sn„) — "--^Go-bulk — W-E'sn-bulk, (1) 

where -Ecci^GcnSn^) is the total energy of the fully re- 
laxed cell with the defect, divided by the total number 
of atoms in the cell. i?Gc-buik andi?sn-buik are the bulk 
Ge and Sn cohesive energies per atom, corresponding to 
perfect diamond lattices of Ge and gray tin (a-Sn). We 
have also calculated the partially relaxed local defect en- 
ergies (E^'^^ ) per cell atom: fixing the cell volume at the 
value Vce corresponding to bulk Ge, allowing only relax- 
ation of the positions of the individual atoms inside the 
cell. 

In order to estimate the elastic field caused by the pres- 
ence of the defect, we also determined the pressure ex- 
erted by each type of local defect on the lattice. This is 
obtained from the dependence of the total energy with 
the volume (i.e. using partially relaxed energies at fixed 
cell volume E^). We define the pressure as: 



Vce dV v=Vg. 

calculated at the unit cell volume which corresponds to 
bulk Ge. In our convention, a positive pressure P will 
mean that expansion of the lattice with respect to bulk 
Ge is favoured. 

The numerical calculations were done using the full- 
potential linearized/ augmented plane wave plus local 
orbital (L/APW-l-lo) method, as implemented in the 
WIEN2K code^i2^i2i. The exchange-correlation effects 
were treated within the GGA (generalized gradient ap- 
proximation) using the Perdew-Burke-Ernzerhof form^^. 
The radii of the muffin tin spheres (Rmt) were chosen 
to be 2.25 Bohr for both Ge and Sn atoms. The cut- 
off parameter RuT-Kmax for limiting the number of the 



plane waves was set equal to 8, with Kmax (= 8/2.25 = 
3.55 Bohr~^) the largest reciprocal lattice vector used in 
the plane wave expansion. To integrate inside the Bril- 
louin zone (BZ) we used a k sampling with (7x7x7), 
(6x6x6) and (4x4x4) Monkhorst-Pack^^ meshes, 
for unit cells with 8, 16 and 64 atoms respectively. In all 
cases, all atoms in the unit cell were allowed to fully relax 
independently and, except otherwise stated, also lattice 
parameter (i.e. cell volume) relaxation was allowed. In 
Table [I] we summarize the results of our ab-initio elec- 
tronic structure calculations. The first column identifies 
each of the local configurations considered. We have la- 
belled the unit cells by explicitly indicating the number 
of atoms of each species they contain, further details will 
be given below. 

First, we have considered the two limiting cases of the 
alloy, namely the perfect diamond lattices of Ge and gray 
tin, in order to extract their cohesive energies as reference 
values for the local defects. In these cases (labelled as 
Ge2 and Sn2 in Table I), we used the Wigner-Seitz cells 
for diamond, containing two atoms and periodic bound- 
ary conditions. We also calculated the corresponding 
bulk energies and lattice parameters for cells including 
64 atoms: Ge64 and S'n64, resulting in cohesive ener- 
gies differing from the 2-atom cell values by less than 0.1 
meV. Observe in the table, that also the lattice parame- 
ter changes are negligible. 

In the 8-atom cells we studied a single substitutional 
a — Sn defect (the GelSnl case of Table I), and the pair 
defect consisting of two contiguous a defects in which 
the two Sn atoms align parallcU to the (111) symmetry 



TABLE L Results of electronic calculations using Wien-2K 
code {L/APW + lo with GGA) for the cells with configura- 
tions described in the text. Successive columns report: corre- 
sponding "alloy concentration" [x); lattice parameter (a) of 
fully relaxed crystal structure; cohesive energy (-Ec) (relative 
to bulk Ge and Sn); cohesive energy per cell atom: fully re- 
laxed local defect energy {Ed) \ pressure (P) due to defect; and 
partially relaxed local defect energy {E^'^") per cell atom: at 
Ge cell volume. All energies given in eV. 
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axis of the diamond lattice (case Ge6Sn2\\). In the latter 
case, we have verified that starting from different orien- 
tations of the pair of Sn atoms, this symmetry is the 
one corresponding to the minimum of the total relax- 
ation process. These 8-atom configurations were studied 
before, with local density approximation for density func- 
tional theory (LDA).— In agreement with the well-known 
general trends, the GGA yields slightly lower cohesion 
energies and larger lattice parameters than LDA, in the 
cases where we can compare our results with previous 
LDA calculations.^^ 




FIG. 1: a — Sn (substitutional) defects in Ge lattice. 
(GelSSnl) 




FIG. 2: P — Sn defect in Ge lattice: 16-site supercell shown. 

{GelASnlp) 



We also studied the single a — Sn defect and the pair 
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FIG. 3: Sn pair-defect in Ge lattice: 16-site supercell shown. 
(0614571211) 



defect configurations in a 16-atom cell, to see the effect 
of the size of the cell on the relaxed configurations. In 
particular, the fact that the energy obtained for two in- 
dependent (non-contiguous) a — Sn defects in a 16-site 
cell {Gel4:Sn2) is approximately twice the energy of one 
a defect in the 8-site cell (GeTSnl) is a further indication 
that our results are reasonably size-independent. 

In Figure[T]we show the the substitutional a—Sn defect 
Gel5Snl, where a 5*71 atom at the center of a 16-site 
cell. In Fig. [21 we depict the non-substitutional (3 — Sn 
defect with a Sn atom at the center of a divacancy in a 
16-site cell of the Ge lattice (denoted GelASnlfj case), 
while G'el4S'n2|| labels the above mentioned pair defect, 
now placed in a 16-site Ge cell (see Fig. All figures 
included in this section, exhibit the atoms in the fully 
relaxed configuration obtained for each cell. 

For 16-site cells we have also studied the effect of ag- 
glomeration of Sn: GellSnb corresponds to five sub- 
stitutional Sn atoms at the center of a 16-site cell; 
G'ellS'n4a denotes the previous configuration with the 
central Sn atom removed, such that a vacancy sur- 
rounded by four Sn atoms occupies the center of the 
16-site Ge-lattice cell, while GellSnAi, refers to the same 
cell containing in its center four Sn atoms plus a vacancy 
now being placed at one of the nearest neighbour sites of 
the cell center. We have also studied two 64-site cells, 
shown in Figs. [3] and [51 containing respectively a 6-5*71 
cluster-defect with a central divacancy (G'e56S'n6), and 
a P defect surrounded by 6 S'7i-atoms {Ge56Sn7[}>), al- 
ways in the Ge lattice. 

The last row of Table I refers to GeSn, the x=0.5 
zincblende structure, which should be strain-free. No- 
tice that the pressure shown in last row of Table I, was 
calculated relative to the bulk-Ge cell volume (Eq. [2]), 
instead of relative to the zincblende cell volume. 

Our ab-initio electronic calculations, performed in su- 
percells larger than in previous studies^i^iSi and with 
various different local defects, among other facts reveals 
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FIG. 4: 6-Sn cluster-defect with a central divacancy in Ge 
lattice: 64-site supercell shown. {Ge56Sn6) 




FIG. 5: 7-Sn cluster-defect (/3 defect surrounded by 6 5*71,- 
atoms) in Ge lattice: 64-site supercell shown. {GebGSnT/ji) 



that: 

i) cubic octahedral symmetry is favoured in clusters 
with few Sn atoms, as can be seen in the fully relaxed 
configurations of Figs. [21 [H and [5l 

ii) electronic energy per atom is not much increased by 
accumulating Sn clusters; 

iii) substitutional a — Sn and non-substitutional 13 — Sn 
single-atom defects do create opposite elastic fields: this 
leads to an effective attraction between them, such that 
if they do merge elastic energy is released and the small 
cluster relaxes to cubic symmetry; 



iv) building up a local pressure in the lattice is not 
necessarily correlated with the electronic energy of the 
cluster. This agrees with former ideas about the exis- 
tence of two independent difficulties for the formation of 
homogeneous alloys, namely, the large size difference be- 
tween Ge and Sn atoms with a 15% lattice mismatch, 
and the d-bands which make gray tin unstable at room 
temperature^Siii^ 

v) the energy of two a — Sn defects increases when 
they are closer: indicating it would be difficult to form 
a homogeneous high concentration totally substitutional 
alloy (due to the electronic energy cost). If, as generally 
mentioned before, strain was the only/main factor jeop- 
ardizing formation: x=0.5 zincblende {GeSn) should be 
easy to form, which is not the case. 



B. Statistical mean field approach 

We used the valuable information from our electronic 
calculations of the previous section, to try to understand 
and predict the behavior of the alloy as a function of the 
Sn concentration (x). Taking advantage of the fact that 
both experiments and calculations indicate a favourable 
substitutional incorporation of Sn in the Ge matrix for 
low concentrations, we could model the statistical be- 
havior of the mixture during annealing to attain ther- 
modynamical equilibrium by a stochastic approach. At 
low x, defects containing only one Sn atom are relevant. 
Imagine a tetrahedral cage of Ge atoms in which the 
central four-coordinated site could be occupied by either 
Ge or Sn, the former representing the perfect Ge lattice 
(Tfje) while the latter represents an a — Sn defect in Ge 
[To). By substitution of only one Sn atom in any of 
the corners of the cage, one gets the 6— fold coordinated 
S'n— pair defect studied in previous section. This 6— fold 
cage, which is the environment of a divacancy, could also 
be occupied by a single central Sn atom, giving the non- 
substitutional /? — Sn defect studied in previous section 
[Tfj). The a — Sn defect can be transformed by substitu- 
tion of a single atom into the pair-defect, which becomes 
important only for larger x. To keep the treatment to 
first order in a;, we approximate the pair defect by two 
separate a — Sn defects, based on the fact that the nu- 
merical calculations of previous section indicate that the 
energy of the pair-defect is quite similar to twice the 
energy of one a — Sn defect, Ea- However, this energy 
difference will be taken into account below, when needed 
in the treatment. 

To first order in x, there are only these three local 
configurations: Ti{i = a,P,Ge), which form a complete 
closed set, in the sense that one can transform one into 
the other without leaving the set by substitution of a 
single atom. This allows one to formulate a dynamical 
process, representing the formation of the infinite solid by 
a large number of such one-step transformations encoded 
in a stochastic matrix (which we will denote by S) .27,28 
The elements of S can be considered as the probabilities 
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Pa,b of obtaining configuration b starting from configu- 
ration a, and must contain all the physical information 
available. Concretely, writing a one-step transformation 
as: 



Svq = Vi, 



(3) 



where vq is a vector whose three components represent 
the initial concentrations of each of the three configura- 
tions Ti, therefore the sum of its three components equals 
one (i.e. normi{vo) = 1). In general, the final vector 
vi is different but its nornii should be conserved. This 
is guaranteed by normalizing the columns of matrix S 
to one, so that S becomes a stochastic matrix. A well- 
known property of stochastic matrices is that one of their 
eigenvalues Ai is equal to 1 while all others have magni- 
tudes which are less than 1 (|Ai| < 1) . After repeating 
the transformation n times, the transformed vector can 
be expressed in terms of the eigenvalues and eigenvectors 
Bi of the stochastic matrix S as. 



We shall consider the remaining transition probabili- 
ties to be zero, since they are second order in y, except 
for the transition between two neighbouring a configura- 
tions described by: 



Pa-c ^ {x- yf 4232 e- 



-(£;||-B„)/(feBT), 



(8) 



due to the energy considerations which were mentioned 
above. Observe that the energy barrier is only the differ- 
ence between the pair defect and the a-defect, since only 
one step is needed to form it from the initial a configu- 
ration. 

Summarizing, the stochastic matrix S takes the follow- 
ing form: 



AO B 
C 

{l~A) 1 



(9) 



where: 



i=l 



(4) 



where ai are the components of Vf) in the eigenvector 
basis. If n — > cxd, the only surviving clement of the 
summation above is: 



S"wo ~ (Ai)"aiei = AiOiei = oiei. 



(5) 



Taking into account that normi{vo) = 1 = normi(ei), 
and that S is a stochastic transformation, we immedi- 
ately obtain: oi = 1. Therefore, the elements of ei 
should be considered as the concentrations of the final 
configurations in the macroscopic system. For any ar- 
bitrary initial concentrations, i.e. for any vq, under ho- 
mogeneous formation conditions, one attains the fixed 
point: 



ei 



(6) 



where y is the concentration of /3— defects in the sample 
(and X — y denotes the concentration of a— defects). 

The transition probabilities Pa,b = P{a,b)/Na,b (Na.b 
are the required normalization factors) are composed of 
three parts. A first factor is the conditional probability 
of obtaining configuration b, starting from configuration 
a, that is: XaXb- The second factor is a number count- 
ing all the possibilities of arranging the chemical bonds 
without changing the final results of the transition, and 
it could be interpreted as a sort of local configurational 
entropy. Finally, the energy barriers to perform the tran- 
sition should be considered. If the time scale for a chem- 
ical bond to reach thermal equilibrium is much smaller 
than the time it takes for the whole solid to equilibrate, 
then the third part can be written as a Boltzmann factor 
containing the temperature of the bath. In our case, we 
thus have: 



Pa- 


Ge 


= (1 


— x){x — 






Pp- 


Ge 


= (1 


— x)y6\e 


-Efil(kBT). 

1 




Pae- 


-Ge 


= (1 


~xf4k- 




(7) 



A 

B 
C 



Ba-Ge 



Ba-Ge + Bp^Ce + Bce-Ge 

Bp-Ge 

Ba-Ge + Bp^Ge + Bce-Ge 



(10) 



Notice that, due to the normalization condition, only 
three independent energies appear in the problem, as can 
be verified by respectively dividing numerator and de- 
nominator of A,B and C by Bce-Ge- Thus, only the 
three relative energies of each defect configuration with 
respect to the Ge lattice (i?^ = i?^ — Eoe) will appear, 
which are precisely the energies of the local defects cal- 
culated numerically in the previous section. 

The eigenvector of S with eigenvalue 1 is: 



ei = {Xa,X/3,XGe), 



(11) 



where: 



B 



1-A 

Xf3 = CxGel 



XGe\ 



I- A 



XGe 



{1-A)il + C) + B' 



(12) 



which due to the homogeneity condition should be iden- 
tical to Eq. [51 

Therefore, at a fixed temperature (T) one can obtain 
the concentration of /3 — Sn defects by finding the zeros 
of the function 



f{x,y) = 



I- A 



{l-A){l + C) + B 



(1-x), 



(13) 



which we do numerically. The needed values of the en- 
ergies are taken from Table U and are: Ea = 0.018 eV, 
Ep = 0.120 eV, and E\\ = 0.041 eV. 
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0.3 0.4 0.5 0.6 0.7 
Sn concentration (x) 

FIG. 6: Concentration of /3 — Sn defects as a function Sn 
concentration in a Gei_a:Sna; alloy, at different temperatures. 
Full line: T = 16°C; 39°C (dotted); 62°C (plus signs), 85°C 
(asterisks), 108°C (crosses); 131°C (circles). 



In Fig. [5] we depict the zeros of f{x,y) obtained for 
six different temperatures. Notice that the only physi- 
cally meaningful solutions are the ones below the y = x 
line (since the /3— defect concentration y cannot be larger 
than the total 5'n— concentration present, x). Observe 
that at room temperature there are no P defects present 
for concentrations x < 0.2, meaning that the alloy is 
perfectly substitutional (only a — Sn defects present). 
Therefore, a continuous shift of the electronic properties 
is to be expected and could be accurately modelled by 
a simple VGA approximation, as previously noted and 
measured^'^" and we will further discuss in next section. 
For larger 5'n-concentration, x > 0.2, one finds that a 
non-zero concentration of /? — Sn defects appear, which 
would introduce electronic traps in the gap of the alloy. A 
monotonous increase of the number of /3 defects follows, 
as more Sn is incorporated to the Ge lattice. 

Furthermore, our simple model predicts that as the 
temperature of formation of the solid is increased the /3 
defects start to appear at lower 5*71 concentrations, as 
expected because the increase in thermal energy favours 
the overcoming of all energy barriers. This fact provides 
an explanation for the findings of Gei-xSn^ alloy prepa- 
ration by chemical vapor deposition (CVD)^, where it 
has been observed that more substitutional Sn could 
be incorporated if the substrate temperature is lowered. 
This is also known to happen in amorphous Gei„a;Sna; 
alloysii. 

In Fig. [7] we show an amplification of the low concen- 
tration results, where the above discussed dependence 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 
Sn concentration (x) 

FIG. 7: Amplification of the low concentration region of 
Fig.El Full line: T = 16°C; 39°C (dotted); 62°C (plus signs), 
85° C (asterisks), 108° C (crosses); 131° C (circles). 
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FIG. 8: Comparison of use of alternative defect energy defi- 
nitions: relaxed vs. fixed at Voe cell volumes, on the concen- 
tration of /3 — Sn defects as a function of total Sn concentra- 
tion. Using Ed, at : T = 18°C (squares); 79° C (asterisks); 
134°C (circles). Using £'J'^% at the same temperatures: 18°C 
(crosses), 79° C (plus signs); 134° C (dots). 



on temperature is clearly seen. Also, notice that for high 
enough temperature the non-physical roots start to cross 
the X = y line, meaning that f3 defects can be formed at 
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any 5'n— concentration. This is expected to happen, be- 
cause disposing of large amounts of thermal energy allows 
the significant creation of more complex defects, than 
those considered in this simple approximation. We think 
that there is not much to gain by increasing the com- 
plexity of the present model, since the larger the defect 
space considered, the less reliable are the local energies 
calculated. 

Finally, in Fig. [8] we show the effect of using the alter- 
native definition for the defect energy obtained by partial 
relaxation, fixing the cell volume at the bulk-Ge value 
(iSj^"). As expected from the differences in respective 
energies stated in Table I, which show a relatively lower 
increase for the /3-defect energy than for the a or pair- 
defect cases, fixing the cell volume leads to /3— defects 
appearing at slightly lower Sn concentration values. The 
comparison of results using the alternative defect energy 
definitions reveals that the main effect of using the partial 
relaxation is a general shift (reduction) of the "effective" 
energy barriers in the Boltzmann factors of Eqs.[7|and [H 
thus being analogous to a rescaling of the temperature. 



C. Electronic structure studies 

1. Supercell calculations 

The ab-initio local defect calculations presented in sec- 
tion nil Al also yield the density of states and energy 
bands of perfectly ordered periodic crystals formed by 
those supercells. As an illustration we show in Figure [5] 
the densities of states (DOS) of three characteristic crys- 
tals, namely the pure Ge and Sn ones, and the zinc- 
blende-type GeSn, previously studied in Refi^. 




one would expect. Notice that the pure Ge DOS does not 
have a gap at the Fermi level, which is not a surprise due 
to inaccuracies at small energies around the Fermi level 
of the approach GGA-Wien2K. Another peculiar thing is 
the fact that deep down the valence band (around -SeV), 
both pure systems present a single Van Hove singular- 
ity with zero DOS corresponding to degeneracy of the 
bands at the X point of the first Brillouin zone of the 
fee diamond structure. This degeneracy is due to the 
inversion symmetry, which is lacking in the zincblende 
structure. The wide gap in the density of states around 
those energies exhibited in zinc-blende GeSn is due to 
the removal of this degeneracy.^ Of course, such fea- 
tures due to symmetry-breaking are not expected in the 
real crystalline alloys, since perfect order is not realized. 
All this can be seen more clearly in Figure [TUl where the 
electronic bands from the 2-atom cell are shown. Observe 
that the lowest conduction band crosses the Fermi level 
at point L. 




X Z W K 



FIG. 10; ab-initio (Wien-2K) calculated bands with 2-atom 
cells for GeSn. 



FIG. 9: ab-initio (Wien-2K) calculated DOS with 16-atom 
supercells. 

It can be noticed in Figure [9] that the Ge bands are 
slightly wider than the Sn ones, and that the density of 
states of GeSn at the Fermi level lies between the higher 
value for Sn and the very small value obtained for Ge, as 



2. TB + VGA 

Tight-binding (TB) calculations for electronic band 
structures are very useful because, besides their simplic- 
ity, it is possible to incorporate many features by suitably 
choosing the TB parameters. These method applied to 
Ge proved not to be straightforward, since a naive cal- 
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culation considering sp^ orbitals is unable to reproduce 
essential features like the indirect band gap and the band- 
widths. This is due to the mixing of d electrons in the 
conduction bands of Ge, which for Si or C are not impor- 
tant, while for Sn and Pb are responsible for their metal- 
lic behavior. The role of d electrons could be mimicked 
by introducing a pseudo-orbital (s*) without increasing 
the size of the TB Hamiltonian. Finer properties, such as 
the exciton spectrum of Ge cannot be explained without 
introducing spin-orbit interactions. 



VCA+TB bands - Ge_(1 -x) Sn x lor x = 0.1 5 




L X U,K 



FIG. 11: TB-hVCA calculated bands as in Ref. ^ for 
Gei_a;Sna; at x=0.15. 

In 1987, Jenkins and Dow^ presented a complete TB 
model for the Gei_j;Sn^ alloy, including spin-orbit, s*, 
and second- neighbour interactions. They examined sev- 
eral properties by using a simple virtual crystal approxi- 
mation (VCA)iS description of the alloy. Their results 
are in very good agreement with experiments for low 
concentrations, in particular they predict the indirect 
to direct gap transition to occur at x « 0.15. For the 
sake of completeness and to facilitate comparison with 
the low concentration regime, we have reproduced their 
TB-)-VCA calculations and in Figure [11] we depict the 
electronic bands precisely at this critical concentration. 
Notice that the optical gap is the same at the center of 
the 1st Brillouin zone and at point Li^ 

The dependence of electronic properties as a function 
of concentration continues smoothly in this approxima- 
tion, contrary to the experimental facts. To illustrate 
this, in Figure [12] we plot the bands at 50% concentra- 
tion, where a small energy gap at the Fermi level is still 
present. Spectroscopic ellipsometry and photoreflectance 
experiments-^ found that there is a strong bowing ten- 
dency, i.e. non-linear behavior of this gap with concen- 
tration, that would close the gap at a much lower value: 
X ~ 0.37. As one expects for the alloy, no degeneracy is 
removed at point X, because it does not introduce spu- 
rious order like the zinc-blende bandstructure of Fig. [10] 

In summary, these TB-I-VCA calculations can be 
trusted only at low concentrations, where they provide 
a valuable tool to interpret experiments. 



VCA+TB bands - Ge x Sn_(1 -x) lor x = 0.5 
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FIG. 12: TB+VCA calculated bands as in Ref. for 
Gei-xSn^; at x=0.5. 



3. Coherent Potential Approximation 

In many random alloys, the coherent potential approx- 
imation (CPA)^ for the treatment of substitutional dis- 
order leads to an improved description of the electronic 
structure, overcoming some limitations of the VGA dis- 
cussed in previous section. CPA has been applied suc- 
cessfully to many semiconductors;^ and in particular to 
Gei_2;Six alloys^ where it predicts moderate bowing in 
the optical transition energies. Concerning Gei-^jSn^:, 
special supercell calculations have been performed to 
mimic the alloy by averaging over selected structures^i^. 
They predict a large and compositional-dependent direct 
gap bowing, in agreement with experimental reports.'^^ 

CPA electronic structure studies of the Gei_a;Sna; al- 
loy were lacking, up to now. This has prompted us to 
undertake them, as a further way to indirectly test our 
hypothesis regarding the relevance of non-substitutional 
defects in this alloy. In general, effective-field calcula- 
tions, like VGA or CPA, give excellent predictions of the 
electronic properties in substitutional alloys, thus a fail- 
ure to describe the experimental findings suggests a fur- 
ther indication of the presence of /3 — 5*71 defects. In fact, 
we found that if one takes into account the structural 
changes undergone as a function of Sn concentration, one 
can extend the range of validity of these approximations. 

In the following, we report results of our CPA calcu- 
lations for the Gei-a^Sn^: alloy, using the FPLO-5 code 
developed by the IFW-Dresden group.^^ The FPLO- 
5 package is a full-potential local-orbital minimum-basis 
code, combining ab-initio LSDA (local spin density ap- 
proximation) bandstructure treatment^^ with single-site 
CPA routines for substitutional disorder— (it employs 
the Blackmann-Esterling-Berk single site-CPA extension 
for combined diagonal and non-diagonal disorder'^^). 

In Fig. [13] we show the FPLO+GPA total den- 
sity of states for Gei-^^Sn^ as a function of energy, 
at different substitutional-S'n concentrations: x = 
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FIG. 13: (Color online) CPA for Gei_^Sn^, via FPLO-5 code: 
total density of states for the a — Sn concentrations indicated 
inside plot. Lattice parameters from respective Wien-2K ob- 
tained structural data. 



125-DOS 




E [eV] 



FIG. 14: (Color online) CPA results for Geo.875Sno.125, via 
FPLO-5: total and partial (species-resolved) DOS for each 
alloy component, as a function of energy. 



0.0, 0.062, 0.125, 0.25, 0.5. Notice that a smooth behavior 
as a function of concentration is obtained, with changes of 
the bandwidth, and a progressive filling of the gap at the 
Fermi level with S'n-concentration, as one would expect. 
Experiments in Gei_a;Sna; alloys^ have indicated that 
the direct energy band gap decreases primarily through 
an increase in alloy concentration (and applied coherency 
strain mainly reduces the valence band DOS, instead of 
the magnitude of the gap). 

The results we present were obtained adjusting the in- 
put lattice parameters at each alloy concentration using 
our respective Wien-2K obtained structural data (Table 
1). If, as it is usual in CPA, one uses the fixed lattice 
parameter corresponding to the perfect (x=0) lattice at 
all concentrations, one finds that the trends with concen- 
tration are much smaller, if observable. For instance, the 



resulting valence band obtained is almost unchanged by 
concentration increases of up to 25%, thus not reflecting 
the effect of the difference of bandwidths (and shapes) of 
the two alloy components. In this sense we have found 
that some improvement is gained by supplementing this 
effective field method with information on the structural 
changes undergone by the alloy as a function of concen- 
tration, considering that we found limitations of the tra- 
ditional CPA description for Gei_a;Sna; even at relatively 
low concentrations. 

To better convey the information obtained by the 
FPLO-I-CPA approach, in Fig. [H we show the CPA re- 
sults at x=0.125, namely the total density of states as 
well as the orbital-resolved partial densities of states. At 
this low Sn concentration, the total DOS follows closely 
the Ge-host one. 

Increasing concentration to x=0.25, in Fig.[Tn]we show 
the total DOS as well as the orbital-resolved partial DOS: 
here the most relevant ones correspond to s orbitals 
(given by the 4s orbitals in Ge, and the 5s ones in Sn) 
and p orbitals {Ap of Ge, and bp of Sn) while, for clarity, 
the relatively smaller d components are not shown. Our 
results show good agreement with the valence band s and 
p partial DOS obtained by a previous LMTO Green func- 
tion study f?- Basically, the spectral weight of p orbitals 
strongly dominates around the Fermi level, while the s 
orbitals are the most relevant at lower energies, for both 
species. 



25- DOS 

25- s 
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FIG. 15: (Color online) CPA results for Geo.75Sno.25, via 
FPLO-5: total density of states, orbital and orbital-fspecies 
resolved main DOS contributions, as functions of energy (ref- 
erences in plot). 

Finally, in Fig. [16] we exhibit an example of the prob- 
lems encountered by CPA, based on the assumption that 
only substitutional a — Sn impurities are present, which 
can not be overcome by taking into account the lattice pa- 
rameter changes as a function of Sn concentration. Here 
the X = \ density of states is shown, in addition to the 
X = 0, 0.5 ones, showing that the smooth trend with Sn- 
concentration of Fig. [T31 consisting of the filling of the 
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E [eV] 

FIG. 16: (Color online) CPA via FPLO-5 calculated DOS for 
Ge, Sn and Geo.sSno.s. 



gap at the Fermi level, is not valid at larger concentra- 
tions {x > 0.5). Contrary to expectations, the DOS at 
the Fermi level for pure Sn is not the largest obtained 
for the alloy series. We believe that this provides a fur- 
ther indication of the presence of non-substitutional Sn 
defects, and their relevance for a complete and consistent 
treatment of the problem. 



IV. CONCLUSION 

We have proposed a mechanism to understand the pe- 
culiar properties of Gei-ajSn^; alloys. The main assump- 
tion is the existence of the p — Sn defect, occupying a 
divacancy in the Ge host, which imposes a severe strain 
in the lattice oppposite to the one caused by the a — Sn 
substitutional defect. The other feature of the (3 defect 
is that it causes a six-fold octahedral coordination on the 
Sn atom, favouring the nucleation of white tin inclusions 
which eventually segregate. 

We demonstrate the plausibility of this assumption by 
performing ab-initio electronic calculations in supercells 
containing different local defects. Starting from the sin- 
gle atom a and (3 defects, and increasing the number of 
Sn atoms, in order to estimate: the electronic energy, the 
relaxed configuration around the defects, and the pres- 
sure caused by them. These studies suggest that the 
cubic octahedral symmetry is favoured in clusters con- 
taining few Sn atoms, and that the electronic energy 
per atom is not very much increased by accumulating 
clusters of Sn. The calculated pressures indicate that 
the large positive elastic field around these large defects 
may be released by the aggregation of f3 defects, which 
have an opposite elastic field. We have observed that 
building up a local pressure in the lattice is not necessar- 
ily correlated with the electronic energy of the cluster. 
This is in agreement with former idea s^°'^^ that there 
are two independent reasons to account for the difficul- 



ties to form homogeneous Gei-a^Sn^, alloys, namely the 
electronic d— bands which make gray a — Sn unstable at 
room temperature and the large difference in the size of 
the atoms. We have also observed that the energy of 
two a defects increases when they are closer together, 
also indicating that it would be very difficult to obtain a 
homogeneous alloy formed only by a defects since their 
mean-distance would be smaller than a unit cell. Notice, 
that if strain was the only factor jeopardizing formation, 
the X = 0.5 zincblende compound should be easy to form, 
since it would be strain-free: but this is not the case. 

The knowledge gained by the local defect studies pro- 
vided the basic ingredients to formulate a simple statis- 
tical model to investigate the relative concentration of 
a and P defects in thermal equilibrium as a function of 
the total Sn concentration, x. The model focuses on the 
dynamical evolution of few local configurations by using 
a stochastic matrix containing all the statistical informa- 
tion to calculate the transition between one configuration 
to another. The results show that, at room temperature, 
there are no /? defects present for x < 0.2. Thus, our 
model supports the experimental finding that, at con- 
centrations below 20%, Sn only enters substitutionally 
with a local tetrahedral environment. We find that the 
number of P defects at a given concentration increases 
with temperature. We also found that the concentration 
at which /3 defects start appearing decreases at larger 
temperatures. This suggests that one could obtain ho- 
mogeneous alloys with higher Sn content if one decreases 
the temperature of the thermal bath. This would mean, 
in the case of epitaxial growing by CVD or other tech- 
niques, that lowering the temperature of the substrate on 
which the alloy is grown favours the formation of homo- 
geneous alloys at higher Sn concentrations. This agrees 
with the experimental findings.^" In particular, we pre- 
dict that a; = 0.5 homogeneous alloys would be possible 
to form at temperatures below -90° C. 

In the last part of our work, we have tried to gain more 
insight into our problem by using effective-field methods 
in order to investigate the low concentration region in 
more detail. There, we found that TB-I-VCA, while pro- 
viding a good description of the substitutional-dominated 
regime, fails to predict the closing of the electronic gap 
at x = 0.37 seen in experiments. Furthermore, more so- 
phisticated mean-field techniques as the CPA (even us- 
ing a combined ab-initio realistic bandstructure with the 
single-site substitutional CPA), does not provide much 
improvement. However, if one supplements the CPA 
by introducing the approximately linear concentration- 
dependence of the lattice parameter, obtained from our 
ab-initio local defect calculations (in agreement with ex- 
periments), a slight improvement of the results is ob- 
tained. All this, further reinforces the hypothesis put 
forward in the present work. We think that an extension 
of the CPA^ able to take into account more complex 
defects, such as the non-substitutional /3 defect, would 
certainly be extremely useful to improve the description 
of the electronic properties of these fascinating and use- 
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ful alloys, and to study their effect on the large bowing 
of the direct band gap and the indirect to direct gap 
crossover Research along these lines is cur- 

rently in progress. 
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